Strict and fussy mode splitting in the tangent space of the Ginzburg-Landau equation 
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In the tangent space of some spatially extended dissipative systems one can observe "physical" 
modes which are highly involved in the dynamics and are decoupled from the remaining set of 
hyperbolically "isolated" degrees of freedom representing strongly decaying perturbations. This 
mode splitting is studied for the Ginzburg-Landau equation at different strength of the spatial 
coupling. We observe that isolated modes coincide with eigenmodes of the homogeneous steady 
state of the system; that there is a local basis where the number of non-zero components of the 
state vector coincides with the number of "physical" modes; that in a system with finite number of 
degrees of freedom the strict mode splitting disappears at finite value of coupling; that above this 
value a fussy mode splitting is observed. 
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INTRODUCTION 

Nonlinear dissipative spatially extended systems have, 
from the formal point of view, infinitely many degrees 
of freedom. But many important examples are known 
where the chaotic solution of an extended system evolves 
in an effective manifold of finite dimension that is called 
the inertial manifold [ij. H. Yang et al. [1| suggest 
that the tangent dynamics of Kuramoto-Sivashinsky and 
Ginzburg-Landau equations is essentially characterized 
by a well-defined set of vectors called "physical" modes 
which are decoupled from the remaining set of hyper- 
bolically "isolated" degrees of freedom. In this case the 
physical modes can be a local linear approximation of 
the inertial manifold, while isolated modes are orthogo- 
nal to this manifold and are responsible only for transient 
processes. 

The structure of the tangent space of a dynamical sys- 
tem is characterized by Lyapunov exponents and associ- 
ated with them Lyapunov vectors. There are two or- 
thogonal sets of vectors called backward and forward 
Lyapunov vectors Q. They can be computed in the 
course of the standard procedure of computation of Lya- 
punov exponents [3| in forward and backward time, re- 
spectively 0, These vectors are not covariant with 
the tangent dynamics in a sense that the tangent map- 
ping, being applied to them, does not produce the for- 
ward or backward vectors. Though the existence of the 
covariant Lyapunov vectors (CLVs) was known for a long 
time, they became available only recently thanks to ef- 
fective numerical algorithms ^ 0] . These vectors are not 
orthogonal, they are invariant under time reversal and 



covariant with the dynamics. Because these vectors span 
Oseledec subspaces, they allow access to hyperbolicity 
properties f6|. CLVs are the generalization of the notion 
of "normal modes". They are reduced to the Floquet 
vectors if the flow is time periodic and to the stationary 
normal modes if the flow is stationary [7| . 

In this paper we study the mode splitting reported in 
Ref. |2|. The motivating idea is very simple. Consider 
an extended dynamical system. When the spatial cou- 
pling is very strong, the effective dynamics should be low- 
dimensional due to synchronization effects. It means that 
the number of physical modes should also be small. But 
when the coupling is very small, the spatial cells become 
almost independent. In this case all degrees of freedom 
are important so that the mode splitting can vanish. We 
study a chain of amplitude equations that appear from 
the Ginzburg-Landau equation when spatial discretiza- 
tion is introduced. The step size of the discretization is 
used as a control parameter. Varying the step we analyze 
the mode splitting at different intensities of coupling. 

The paper is organized as follows. In Sec. U the model 
system is described. Sec. |ll] is devoted to the case of 
strong coupling when the system is close to the bifurca- 
tion point. In Sec. lIIII the strict mode splitting is analyzed 
that is observed at moderate values of coupling. Sec. IIVI 
represents the case of a weak coupling when the strict 
splitting disappears and a fussy splitting is observed in- 
stead. Finally, in Sec. [Vlwe summarize the obtained re- 
sults. 



THE MODEL SYSTEM 
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Consider a ID complex Ginzburg-Landau equation 
dta — a — {1 + ic)|apa + (1 + ib)d^a. To find solutions of 
this equation numerically, we represent the second spatial 
derivative as a finite difference. In this way, the partial 
differential equation is transformed into a chain of N am- 
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plitude equations: 

a-n = an - (1 +ic)|a„pa„ + (1 + ib)n{an) /h"^ , (1) 

where a„ = a„(t) (n = 0, 1, . . . , — 1) are complex vari- 
ables, /i is a step size of the discretization, and c and 
h are real control parameters. We set c = 3, 6 = — 2 
which corresponds to the regime of so called " amplitude 
turbulence" Q. Function K(a„) determines the diffu- 
sive coupling and no-flux boundary conditions: K(a„) = 
a„_i-2a„+a„+i {n = f , 2, . . . A^-2), «;(ao) = 2(ai-ao), 
K{aM-i) — 2(ajv-2 — ojv-i)- We are interested in the 
properties of this system at different strength of spatial 
coupling. So, the step size h shall be our control param- 
eter. Treating the discrete space representation of the 
Ginzburg-Landau equation as a chain of oscillators al- 
lows us to freely change the step size h without taking 
care of the validity of the numerical scheme. 

To understand what happens when h tends to zero, 
we perform the rescaling in Eq. ^ h eh' , a — >■ a/e, 
and t e^t. In the resulting equation d„ ~ e^a„ — (1 + 
ic)|a„pa„ + (1 + i&)K(a„)/(/i')^ the decreasing of e cor- 
responds to the decreasing of h in Eq. ([T]). The e here 
can be treated as a bifurcation parameter, controlling the 
stability of the homogeneous steady state Q . This state 
becomes unstable at e = 0, and the system enters the 
regime of spatio-temporal chaos at e > 0. So, returning 
to Eq. ([T|), we can say that when h is small the system 
is just a little bit above the point of the emergence of 
spatio-temporal chaos, and it has only a few positive 
Lyapunov exponent. Increasing h results in chaotic dy- 
namics with an increasing number of positive Lyapunov 
exponents. 

II. ISOLATED MODES AND EIGENMODES 

Consider covariant Lyapunov vectors £i of the sys- 
tem ([T]) . When h is decreased and the system approaches 
from above the bifurcation point where the homogeneous 
steady state becomes unstable, CLVs converge to eigen- 
modes of this homogeneous steady state (in fact, these 
are the modes of Fourier decomposition of the solution). 
For no-flux boundary conditions the eigenmodes read 
gmin) = Sm(i)7m(i) cos(fcimn), where ki = tt/{N - 1), 
and m = 0, ±1, • • • ± (A^ - 2), iV - 1. The total num- 
ber is 2N — 2, but because cosine is an even function, 
modes m and — m are identical and only N modes with 
TO > can be considered. Sm{t) is a normalizing factor: 
X)n3™('*)^ — 1- lm{t) is a vector having two compo- 
nents which can be computed using ti. Vector ti has N 
elements {£i)2n corresponding to Re(a) and also N ele- 
ments {ii)2n+i for Im(a). At the bifurcation point each 
£i coincides with one of the eigenmodes, say to. It means 
that dividing {£i)2n and {£i)2n+i by cos(fciTOn) we ob- 
tain a set of N identical couples that are the components 
of 7m (i)- Beyond the bifurcation point, these couples 
are not identical, and we define the vector 7m (0 as the 
average of them. 
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Figure L (color online) Average cosines of the angles between 
CLVs £i (horizontal axis) and eigenmodes (vertical axis), 
(a) h = 0.1; (b) /i = 0.5, i = 42 for the vertical dashed line 
and m = 21 for the horizontal one; (c) h — 0.8, i = 77, 
m — 33. 



If the system is not far from the bifurcation point, £i 
should not diverge too much from g^- To verify this, we 
compute scalar products of each £i with each for many 
time steps and find the average values. Both £i and gm 
are normalized, hence, the scalar products are equal to 
cosines of the angle between corresponding vectors. Two 
vectors of unit length coincide when the cosine is equal 
to 1. 

Figure [T]^ a) show average cosines at ft, = 0.1 when the 
system is close to the bifurcation point and has only one 
positive Lyapunov exponent Ai « 0.084. A large part 
of the diagram is occupied by the black points along the 
diagonal surrounded by white area. (In fact, there are 
pale squares off the diagonal, but this is a numerical ar- 
tifact.) It means that corresponding £i indeed coincide 
with gm- Notice that the points are grouped pairwise. 
This is a manifestation of the above mentioned degener- 
acy of eigenmodes with to and —to. The non-degenerated 
modes are orthogonal to each other and are referred to 
as isolated modes [1]. 

The degeneracy of isolated modes is associated with 
the degeneracy of eigen modes that, in turn, depends 
on the geometry of the system. In our case the no- 
flux boundary conditions leads for the Ginzburg-Landau 
equation to a degeneracy of the order two, while in 
periodic boundary conditions give rise to a four-fold de- 
generacy. 

There is an area in the left bottom corner of Fig. [TJ^a) 
where CLVs differ significantly from the eigenmodes. In 
Ref. [2|] such kind of vectors has been called "physical" . 
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We shall refer to them as active vectors. These vectors 
are discussed in the following section. 



III. STRICT MODE SPLITTING 
A. Angles with eigenmodes 

The number of active vectors at /i = 0.1 is small be- 
cause the system is close to the bifurcation point. In 
Fig. [Ijb) h = 0.5, and the system has 9 positive Lya- 
punov exponents. We observe now a large area of active 
vectors that is clearly separated from the set of isolated 
vectors. The isolated vectors are represented by the diag- 
onal structure. The diagonal is not so sharp as in panel 
(a), which means that now angles between isolated vec- 
tors and corresponding eigenmodes, though small, are 
not equal to zero. Correspondingly, these vectors are not 
quite orthogonal to all other eigenmodes. But neverthe- 
less, the isolated vectors remain very close to the eigen 
modes. The split between isolated and active modes is 
marked by the vertical dashed line at z = 42. Also, the 
area of active modes is bounded from above: the hori- 
zontal dashed line is drawn at rn = 21. It means that 
the active vectors have relatively small angles only with 
"their own" eigenmodes, i.e., with eigenmodes with num- 
bers corresponding to the active vectors. The angles with 
the other eigenmodes are much higher. Thus, the set of 
active vectors span approximately the same subspace as 
the corresponding amount of the eigenmodes. 

There is another non-trivial structure at the right top 
corner of Fig. [ijb). The nature of this area is unclear 
yet, but we conjecture that it consists of active vectors 
that becomes relevant when time is reversed. 



B. Fraction of DOS violation 

The isolated modes do not have tangencies with the 
active modes f2l| . The method of detection of this strict 
mode splitting, suggested in Ref. (3, employs a concept 
of domination of Oseledec splitting (DOS) [9|,[l3l. We re- 
call that for almost every time every vector in the tangent 
space Si(t) of a dynamical system grows asymptotically 
at rate given by the first Lyapunov exponent Ai except 
those belonging to a set S2(t) of measure zero. Similarly, 
almost every vector in 52 (t) asymptotically grows at rate 
A2 except those belonging to a set 53 (t) of measure zero 
relative to 52 (i), and so on. Collection of sets Si{t) em- 
bedded one into another is called the Oseledec splitting 
of the tangent space. The splitting is called dominated if 
each Oseledec subspace is more expanded than the next, 
by a definite uniform factor. Let Xi{t, r) be the i-th local 
Lyapunov exponent, computed at time t and averaged 
over an interval r. The Oseledec splitting is dominated 
at i if Ai(t, r) > Xj{t,T) holds for all i > i, and for aU 
t with r larger then some finite tq [3]. In particular. 
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Figure 2. (a) Fraction of DOS violation at r = 51.2 vs. num- 
ber of Lyapunov exponent, (b) Average projection of the 
state vector onto vector-columns of Q vs. the number of vec- 
tor. Dotted vertical lines marks, from left to right: index of 
the first of two zero Lyapunov exponents i = 10, Kaplan- 
Yorke (Lyapunov) dimension rounded up to the next integer 
i = 19, and the number of active vectors i = 42. 



domination implies that the angles between the Oseledec 
subspaces are bounded from zero [lO]. 

Employing the ideas of numerical verification of DOS 
suggested in Refs. we define the fraction of 

DOS violation in the following manner. Fix an inter- 
val T and compute Xi(t,T) for some time t. The vi- 
olation of DOS takes place if Ai(i, r) < Xj{t,T) for 
J > i. Thus, for each i we check this inequality at 
j > i and add 1 to the i-th site of an array if it holds at 
least ones. Repeating this procedure for different times 
and performing an averaging we obtain the fraction of 
DOS violation at r, which formally can be defined as 
i^r^ — (0 (maxj>i [Aj(t, r) — Ai(t, t)]))j, where Q{z) is 
the step function and < . . . >t denotes the time average. 

In the case of multiplicity A; ~ A^+i, corresponding 
i^r*'' is close to 0.5 because of fluctuations of local Lya- 
punov exponents due to a numerical noise. Otherwise 
Ut decays to zero as r grows. In general, the decay is 
asymptotic, but if the splitting is dominated at i, the 
corresponding fraction z/r*' vanishes at finite tq. Unfor- 
tunately, there is no a well-grounded algorithm of compu- 
tation of To except the straightforward observation of 
as a function of r. Because the law of the decay of v)- 
is unknown a priori, there is no idea how to extrapolate 
of Vt"^ to zero to verify if a finite tq exists. But, anyway, 
in points of dominated splitting v)- decays much faster 
against r than elsewhere. It means that if r is sufficiently 
large, a graph v)- against i provides relevant information 
about locations of the splitting. 

Figure [Ha) shows the fraction of DOS violation at 
h = 0.5. The sharp minimum of v^'^ at i — 42 coin- 
cides with the position of the splitting found in Fig. [Tib). 
This minimum, presuming the vanish of at a finite r, 
means that the active modes are hyperbolically isolated 
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Figure 3. Lyapunov spectrum at h — 0.5. The horizontal 
dashed line mark = 0, and the vertical ones are drawn at 
the same positions as in Fig. [2] 

from all the rest ones. The active vectors, located to the 
left from the splitting point, have sufficiently high Vr"^ . In 

(i) . 

this case vanishes only asymptotically which, in turn, 
indicates frequent tangencies between the active vectors. 
The isolated vectors are represented by a series of sharp 
minima and maxima with the period 2. Above we have 
shown that these vectors at ft = 0.5 are very close to the 
eigenmodes gm- Because of the degeneracy, the modes gm 
and g-m have identical growth rates. Hence, there are 
couples of corresponding isolated vectors £i with identical 
growth rates. In turn, this implies the multiplicity of cor- 
responding Lyapunov exponents. Indeed, the spectrum 
of Lyapunov exponents demonstrates a stepwise behavior 
to the right from the splitting point j = 42 and the step 
is 2, see Fig. |3l (The stepwise structure of the Lyapunov 
spectrum corresponding to isolated vectors was also re- 
ported in Ref. 12] •) Thus, the maxima of Vt'^ in Fig.[5l^a) 
are associated with this multiplicity. The deep minima 
indicate the absence of tangencies between isolated vec- 
tors because of the orthogonality of corresponding eigen- 
modes. 

The curve of i^i*^ in Fig. [2][a) demonstrates two more 
interesting features. We observe another point of split- 
ting at i = 19 where Vr has a very deep minimum. 
The Kaplan- Yorke (Lyapunov) dimension in this case is 
Dl ~ 18.2. Thus, we can conjecture, that there are 
two types of active modes that are hyperbolically iso- 
lated from each other, and the amount of the first type 
modes is equal to the Kaplan- Yorke dimension rounded 
up to the next integer. 

Also notice the sharp spike at i = 10. To explain the 
emergence of this spike we need to recall that the complex 
Ginzburg-Landau equation with no-flux boundary con- 
ditions has two zero Lyapunov exponents [14| . There are 
9 positive exponents and Aio = An = 0, so the spike at 
i = 10 indicates the multiplicity of two corresponding £i. 

C. Projections of the state vectors 

The split of CLVs onto active and isolated vectors 
should have a clear and visible manifestation in the dy- 
namics of a system. The idea that the number of ac- 
tive vectors i is an effective dimension of the system. 



suggested in 0, presumes that there exists a local ba- 
sis where the state vector has only i nonzero compo- 
nents. We consider projection of the state vector of the 
system at time t on the local basis composed of back- 
ward Lyapunov vectors. These vectors are orthogonal 
to each other and span the same subspaces as CLVs . 
The backward Lyapunov vectors can be computed much 
faster then CLVs in the course of the standard proce- 
dure of computation of the Lyapunov exponents [4|, be- 
cause the columns of orthogonal matrices Q converge to 
them |3, 55 . So, we multiply transposed matrices Q"^ 
by corresponding state vectors, accumulate absolute val- 
ues of obtained projections, and then average them over 
large number of steps. (In general case, the homogeneous 
steady state should be subtracted from the state vector 
before the multiplication to avoid an unnecessary shift. 
But for our system the homogeneous steady state is 0.) 

The average projection denoted as sq is shown in 
Fig. mb). We observe the curve that agrees good with 

Vr . Projections onto active vectors are large, while they 
are almost zero for the isolated vectors. These two parts 
of the curve are clearly separated exactly at z = 42. This 
clarifies the nature of the discussed mode splitting: We 
indeed observe that the number of active vectors can be 
an effective dimension of the system. Also notice that the 
largest component of the projection has an index which 
is equal to the Kaplan- Yorke dimension rounded up to 
the next integer, as marked by the middle dashed line in 
the figure. 

The idea of using backward Lyapunov vectors in di- 
mension reduction methods has already been considered 
and rejected as non-promising in Refs. [isl - fTTj . Indeed, 
the validity of this basis is not so obvious. Reasoning 
formally, one can imagine an attractor that has inappro- 
priate orientation in the phase space so that the decom- 
position fails to give correct result. But on the other 
hand, let us assume that we have a small spherical cloud 
of points that surrounds a homogeneous steady state in 
the phase space. When the points evolve, the cloud is ex- 
tended along the most unstable manifold. In this case the 
most information about dynamics is positively carried by 
the first several CLVs. We can guess, that this property 
survives later producing the split into active and isolated 
vectors. Anyway, at least for the Ginzburg-Landau equa- 
tion this decomposition gives very appropriate informa- 
tion concerning the mode splitting. 

The picture illustrated in this section for h — 0.5 is 
quite generic. We can observe the strict splitting of active 
and isolated modes as well as the Kaplan- Yorke mode 
splitting for a wide range of h. But when h becomes 
sufficiently high, the situation becomes quite different, 
as will be discussed in the following section. 



IV. FUSSY MODE SPLITTING 

As discussed above, exactly at the bifurcation point 
CLVs coincide with eigenmodes and all of them are iso- 
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Figure 4. Number of active vectors against the coupling pa- 
rameter h. 



lated, while the set of active vectors is empty. There 
are 2N — 2 isolated modes. When h grows, the isolated 
modes are converted into active ones and this conversion 
occurs at both ends of the spectrum. (Compare small 
structures at the ends of the main diagonal in Fig. [Tfa) 
and the large areas in Fig.[TJb).) At the left end the iso- 
lated modes contribute to the set of active modes, while 
at the right end they fill up the other set of modes, which 
are, conceivably, relevant when time is reversed. Figure|3] 
shows that the number of active modes depends on h al- 
most linearly. Because we have a finite number of modes, 
there is a finite h for which all isolated modes are con- 
verted so that the splitting vanishes. If the conversion 
at both ends of the spectrum takes place symmetrically, 
then the isolated modes disappear when there are — 1 
active modes. In Fig. 2] the splitting of modes indeed 
disappears when the number of active modes is close to 
64 at iV = 65. 

The value of h where the mode splitting disappears 
depends on the number of eigenmodes, that, in turn de- 
pends on the number of oscillators in the chain N. Tak- 
ing into account almost linear dependance of the number 
of active modes against h we conclude that in a chain 
with an infinite number of oscillators the mode splitting 
vanishes at infinite h. Thus, the continuous system can 
have the mode splitting at any strength of the spatial 
coupling. In particular, it gives a criterion of correctness 
of a chain as a model of continuous system: the chain can 
model a continuous system only if the step size is below 
the point where the mode splitting vanishes. 

Now we consider the tangent space above the point 
of mode splitting vanishing. Figure [Ijc) demonstrates 
angles between CLVs and eigenmodes a.t h — 0.8 (there 
are 17 positive Lyapunov exponents in this case). One 
can see that this figure differs much from the panels (a) 
and (b). Sets of vectors are still distinguishable, but their 
boundaries are not strict. The isolated vectors are absent 
at all. The vertical dashed line marks an approximate 
boundary between two clusters of vectors, and the hori- 
zontal one separate the spectrum of eigenmodes onto two 
halves. The vectors from the left cluster have relatively 
small angles with the first half of the spectrum of eigen- 
modes, while the right cluster contains the vectors that 
have relatively small angles with the second half of the 
spectrum. In particular, it means that tangencies within 
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Figure 5. Same as Fig. [2]at h — 0.8. The vertical dashed lines 
are plotted at i = 17, i = 30 and i — 77. 

the clusters occur more often than between the represen- 
tatives of two clusters. It can be treated as a fussy mode 
splitting: The left cluster is preliminary relevant to the 
forward time dynamics, while the right one, though also 
involved, but does not include much. In the reverse time 
the roles of the clusters are exchanged. 

Figure [5] shows the fraction of DOS violation and av- 
erage projections of the state vector onto backward Lya- 
punov vectors at h = 0.8. There is no series of sharp 
minima and maxima representing isolated modes, and 
also there is no sharp step in the curve of projections. 
But the curve of projections indicates that the state vec- 
tor still has preferable directions inside the left cluster of 
CLVs, whose boundary is marked by the vertical dashed 
line. In the other words, the number of vectors in the left 
cluster can be considered as approximate effective dimen- 
sion of the system. The projections to vectors from the 
right cluster is much smaller. 

Notice also that two other features of these curves sur- 
vive. We still can see the minimum of Ur"^ and the max- 
imum of the projection corresponding to the Kaplan- 
Yorke dimension rounded up to the next integer, as well 
as the spike, indicating degeneracy associated with two 
zero Lyapunov exponents. 

V. SUMMARY 

We studied the splitting of modes of perturbations, 
represented by the covariant Lyapunov vectors, into sets 
of active and isolated modes 0] ■ We considered a chain 
of amplitude equations obtained from the Ginzburg- 
Landau equation by substitution of the second spatial 
derivative with its finite-difference representation. The 
size of the step of spatial discretization was used as a 
control parameter while the number of oscillators in the 
chain was held constant. When the step size is asymp- 
totically small, the system approaches from above the 
bifurcation point where the homogeneous steady state 
looses its stability, while increasing of the step results in 
more independent dynamics of oscillators. 
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At the bifurcation point there are no active modes. 
All modes are isolated and coincide with eigenmodes of 
the homogeneous steady state. Their spatial structure 
is determined by the number of oscillators and boundary 
conditions. When the system leaves the bifurcation point 
as the step size grows, the isolated modes are converted 
into active ones so that the number of active modes grows 
linearly with the step size. 

For the considered system, the backward Lyapunov 
vectors was shown to be an appropriate basis where a 
number of essential components of the state vector is 
equal to the number of active vectors. In other words, 
the number of active vectors indeed plays the role of an 
effective dimension of the system, as conjectured in [2| ■ 

The active modes were found to be split into two sub- 
sets that are hyperbolically isolated from each other. The 
coordinate of the splitting point is equal to the Kaplan- 
Yorke dimension rounded up to the next integer. We 
conjecture that this indicates the existence of two types 
of active modes. The nature of these different types is 
unclear yet, and more studies are required. 

At the right end of the spectrum we observed another 
set of modes which is similar to the set of active modes. 
Its nature is unclear yet, but we conjecture that these 
modes become relevant when time is reversed. 

At a certain finite value of the step size the strict mode 
splitting disappears. Because this value depends on the 
number of oscillators in the chain, the vanish of the split- 



ting occurs only for a system with finite number of de- 
grees of freedom and probably can not be observed, in 
particular, for continuous systems. It can be used as an 
estimation of the maximum step size of the spatial dis- 
cretization. If a continuous system has the strict mode 
splitting and its discrete model does not have it, it means 
that the step size is too large. 

Above the point where the splitting vanishes the 
spectrum of modes contains two clusters without strict 
boundaries. This can be treated as a fussy mode split- 
ting. The first cluster contains formerly active modes, 
while the other one corresponds, apparently, to the 
modes mainly involved when time is reversed. The pro- 
jection of the state vector on the backward Lyapunov 
vectors indicates that the number of vectors in the first 
cluster could be an approximation of an effective dimen- 
sion of the system. 
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